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EQUATION  OF  MOTION  OF  AN  AIRCRAFT  WITH  A  JET  ENGINE 

By 

G.  M.  Trakhtenberg 
Introduction 

Contemporary  aircraft  building  is  characterized  by  an  even  larger 
application  of  automation  for  flight  control.  In  the  method  of 
established  rectilinear  flight,  in  the  method  of  takeoff  and  landing 
the  automatic  pilot  is  widely  used.  Calculation  of  the  characteristics 
of  the  automatic  pilot  is  done  on  the  basis  of  a  contemporary • solution 
of  equations  of  the  dynamics  of  the  aircraft-automatic  pilot  control 
system. 

For  the  collective  solution  of  equations  of  the  dynamics  of  an 
aircraft  and  autopilot  equations  are  usually  used  from  the  motion  of 
a  "solid"  aircraft  in  a  mobile  system  of  coordinates.  However  for 
aircraft,  the  mass  of  which  greatly  changes  (as,  for  example,  aircraft, 
supplied  with  liquid-fuel  rocket  engines,  solid  fuel  Jet  engines) , 
use  of  equations  of  a  "solid"  aircraft  for  analysis  of  the  dynamics 
of  a  regulated  system  is  a  rough  approximation. 

Ignorance  of  the  influence  of  variability  of  mass  on  the 
dynamics  of  the  aircraft  sometimes  leads  to  errors  in  the  engine 
arrangement  on  the  aircraft.  Variability  of  mass,  as  this  is  seen 
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from  this  article,  can  render  significant  influence  on  the  damping 
characteristics  of  the  aircraft.  And  besides,  these  characteristics 
depend  on  the  location  of  the  variable  masses  and  the  Jet  engine 
to  the  center  of  gravity  of  the  aircraft. 

Numerical  calculations  show  that  the  damping  influence  of 
variable  masses  in  reference  to  its  own  damping  of  the  aircraft 
can  reach  separate  sections  of  the  flight  trajectory  up  to  300$.  Thu3, 
for  example,  during  take  off  of  the  aircraft  with  the  help  of  an 
Accelerator  its  own  damping  coefficient  with  respect  to  axis  yt  is 
equal  to  nyy  =  -0.014,  and  the  additional  damping  coefficient  from 
the  accelerator  and  turbojet  engine  is  equal  to  Am^y  =  -0.032. 

Usually  at  the  Joint  solution  of  equations  of  a  closed  aircraft-auto¬ 
pilot  control  system  we  do  not  consider  additional  damping  from 
variability  of  masses.  This  circumstance  can  lead  to  errors  in  the 
selection  of  transfer  coefficients  of  the  automatic  pilot  at  separate 
conditions  of  aircraft  flight. 

Listed  in  this  article,  the  equations  of  aircraft  motion  are 
universal  and  applicable  for  aircraft  of  variable  and  constant  mass. 

Placement  of  the  problem 

Jet  engines  create  large  thrust  at  the  expense  of  ejection  of 
large  masses  of  fuel  mixture  in  a  unit  of  time.  Consequently,  mass 
and  moment  of  inertia  of  the  aircraft  significantly  change  in  time. 
With  these  circumstances  it  is  impossible  to  disregard  in  the  equa¬ 
tions  of  movement  of  the  aircraft  the  variability  of  mass  and  moment 
of  inertia. 

For  a  conclusion  of  the  equation  of  motion  of  the  aircraft  it 
would  have  been  possible  to  use  by  equations  of  motion  of  a  rocket 
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[l].  However,  due  to  the  peculiarities  of  the  approach  to  the  solution 
of  the  placed  problem  and  in  connection  with  the  newly  introduced 
conceptions,  we  will  list  the  main  hypotheses  and  computations  of  the 
conclusion  of  the  equations  of  motion  of  the  aircraft  of  variable 
mass. 

The  most  general  case  of  dynamics  of  an  aircraft  of  variable 
mass  is  the  motion  of  an  aircraft  with  turbo-jet  engines  ( TJE) ,  in¬ 
asmuch  as  side  by  side  with  the  separation  of  gases  from  exhaust 
nozzles  takes  place  continuous  connection  of  air  masses. 

However,  not  losing  community  of  proof,  let  us  consider  in 
the  beginning  the  dynamics  of  a  more  simple  case  -  an  aircraft  with 
a  liquid-fuel  rocket  engine,  but  then  we  will  spread  the  received 
result  also  in  the  case  of  an  aircraft  with  a  turbo-jet  engine. 

A  liquid- fuel  rocket  engine  works  on  a  fuel  mixture,  which 
consists  of  a  liquid  fuel  and  a  liquid  oxidizer,  located  in  special 
tanks,  connected  to  the  liquid-fuel  rocket  engine.  Mass,  moment  of 
inertia  and  center  of  gravity  of  the  aircraft  change  at  the  expense 
of  consumption  of  the  fuel  mixture.  In  the  process  of  flight  the  mix¬ 
ture  proceeds  in  the  engine  and  after  combustion  is  ejected  from  the 
engine  with  great  speed  which  also  stipulates  the  appearance  of  reac¬ 
tive  force.  This  reactive  force  is  determined,  practically,  by  the 
taking  of  experimental  characteristics  on  a  testing  stand  In  an 
established  method.  However  work  of  the  engine  on  the  aircraft 
together  with  the  entire  system  of  supplying  of  the  fuel  mixture  is 
distinguished  from  the  conditions,  with  which  these  characteristics 
are  taken  on  the  stand. 

The  engine  with  the  whole  system  of  supplying  of  the  fuel  mixture 
participates  in  all  evolution-  of  the  aircraft  and  its  characteristics 
in  dynamics,  certainly,  more  general,  than  the  characteristics,  taken 
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Isolated  from  the  aircraft  in  statics 


Fig.  1.  Diagram  of  a  liquid-fuel  rocket 
engine:  a)  tanks  with  fuel  mixture,  b) 

wiring  from  tanks  to  engine,  c)  engine. 

In  figure  1  a  liquid-fuel  rocket  engine  is  shown.  If  we  revolve 
the  aircraft  around  axis  Yu  then  with  this  in  the  engine  and  system 
of  supplying  of  the  fuel  mixture  appear  transverse  force  and  moment 
due  to  angular  velocity.  Therefore,  in  addition  to  the  characteristics 
of  the  engine,  taken  on  the  testing  stand,  we  will  consider  the  forces, 
appearing  in  the  system  of  supplying  of  the  fuel  mixture  of  the  air¬ 
craft  during  its  rotation. 

The  main  question  which  is  the  subject  of  our  Investigation, 
is  the  question  of  calculation  of  the  dynamics  of  overflowing  of  the 
fuel  mixture  on  the  aircraft  (  including  the  engine)  and  the  variability 
of  its  mass  on  the  motion  of  the  aircraft  on  the  whole.  For  solution 
of  this  problem  we  will  Introduce  some  idealization,  which  would 
reflect  the  main  problem,  disregarding  secondary  Insignificant  ;henomena . 

The  mass  of  the  aircraft  consists  of  two  parts: 

the  first  -  constant  and  the  second  -  variable  ( the  fuel  mixture) , 
changing  in  flight  with  time  (Fig.  l) . 

In  reference  to  the  constant  mass  the  known  motion  equations  are 
applicable  in  the  system  of  coordinates  connected  with  the  aircraft. 


Fig.  2.  Idealized  diagram  of  a  liquid- 
fuel  rocket  engine.  l) fuel;  2)  oxidizer. 

In  relation  to  the  second  part  (variable  mass)  it  is  necessary 
to  preliminarily  derive  the  equations  of  motion.  In  the  conclusion 
of  the  equations  of  motion  of  the  aircraft  it  is  assumed  that  the  law 
of  motion  of  the  fuel  mixture  in  the  combined  system  of  coordinates 
is  known  or  given  beforehand. 

Taking  the  fuel  mixture  as  a  "solid  mass"  into  every  moment 
of  time,  we  will  compose  the  following  idealized  calculated  model  of 
variable  mass  of  an  aircraft  (Fig.  2)  .  The  variable  mass  (fuel  mix¬ 
ture)  consists  of  "solid  bricks",  moving  from  the  tanks  to  the  engine. 
In  the  engine  every  "brick",  after  the  process  of  combustion,  is 
separated  from  the  aircraft  with  great  speed,  in  consequence  of  which 
appears  a  reactive  force.  With  respect  to  all  remaining  "bricks" 
of  the  variable  mass  of  the  aircraft  are  applicable,  obviously, 
equations  of  the  dynamics  of  a  point,  taking  into  account  the  circum¬ 
stance  that  the  "bricks"  in  the  combined  system  of  coordinates  have 
relative  velocities  and  accelerations. 

Consequently,  the  constant  mass  of  the  aircraft  can  be  considered 
also  as  variable,  the  elements  of  which  have  zero  values  of  relative 
velocity.  Therefore  in  the  future  we  will  not  divide  an  aircraft 
into  constant  and  variable  parts,  but  will  consider  it  on  the  whole 
as  a  variable  mass. 
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Dynamics  of  an  aircraft  of  variable  mass  with  a  liquld- 


fuel  rocket  engine 


On  the  basis  of  the  idealized  model  built  by  us  we  will  consider 
the  mass  of  an  aircraft  gradually  decreasing,  consisting  of  separate 
"solid"  elements,  having  relative  velocities  and  accelerations  in  a 
combined  system  of  coordinates. 

To  the  aircraft  is  applied  a  thrust,  equal  to  the  reactive  force, 
determined  experimentally  on  the  testing  stand. 

Taking  these  positions  as  the  base  of  our  further  reasonings 
we  will  drive  equations  of  motion  of  an  aircraft  of  variable  mass. 

For  a  conclusion  of  general  equations  of  motion  of  an  aircraft  we  ex¬ 
pediently  use  the  combined  system  of  coordinates,  adopted  in  aero¬ 
dynamics. 

Let  A  be  any  fixed  point  of  the  aircraft,  and  B  another  point 
on  the  aircraft,  moving  relative  to  point  A  (Fig.  3). 

We  will  take  an  immobile  system  of  coordinates  Ox  y  z  ,  and  a 

e  o  u 

system  of  axes,  connected  with  the  aircraft  Axjy^i.  We  will  unite 
with  straight  points  0,  A  and  B.  Having  introduced  0A,  OB  and  AB 
corresponding  to  vectors  p  ,  p  and  r,  we  will  write  the  following 
equality:  p  =  pa  +  r.  Taking  the  geometric  time  derivative  from 
vector  p,  we  will  obtain  the  vector  of  speed  of  point  B: 

V  ^  -i- (  l) 

at  dt  ill 


Expressing  vector  V  by  a  projecti  m  on  the  connected  axes  of  coordi¬ 
nates,  during  rotation  of  the  air -raft  with  angular  velocity  o>,  we 


will  obtain  [3]: 


Vx,  —  Vajl,  +  2,  -  y,  -t  ~  . 


Vyt  ■=  V Ay,  “T  UJ.-I  —  ' 


at 


£> 


where 


V,,  =»  l f At,  +  v,  —  A'i  -r  -jj- .  ( 2) 

l/Al  ,  Va,,,  Vm,  —  are  projections  of  absolute  velocity  of  point 
A  on  the  connected  axes  of  the  aircraft. 


„  >»v, ,  >»,, _  are  projections  of  angular  velocity  of  the 

aircraft  on  the  connected  axes  of  the  aircraft, 
‘ill  liii.. _  are  projections  of  relative  velocity  of  point 

ill  '  ill 

B  on  the  connected  axes  of  the  aircraft. 


dx,  d\\  ils  1 
ill  ‘  ill  ’  ill 


ifXY. 


X  / 

'  ✓  cJ 


/  /I  X- 


— *, 


Fig.  3.  Ground  and  connected  systems  of 
coordinates  of  an  aircraft. 


Having  designed  the  vector  of  acceleration  w  on  the  connected 
axes  of  coordinates,  we  will  find  [3] 5 


dl  dt 


dui . 

z, - -£-y,  t  2» 


-2..,  ^-  +  £^4- 

1  <//  dP 


-r  •r  »/,  <■>.«.  l’|  —  *?,  —  *"/,  Vay,  —  <"I,  xi  +  “>/,  “r,  zi  - 

<1Va,,  ,  d-t,  .  rf-T.  .  .  „  </A,  c  d: ,  <Py, 

*  rf/  dl  '  dt  '  dl  dt  dt 1 

+  ‘•j,  Vaj,  +  id,,  .r,  —  y2  —  id,,  Vaz,  —  xij.y,  +  (dx, 

rf*.  rf«r,  .  ,  „  i/y,  „  ,/.t,  ,  d-z,  , 

—  —  -f - y.  —  —  -  a.  t  2wj.  —  —  2<*v. - h  " — 

dt  dt  dt  ‘  '  dt  'lit  dO 

-r  “a,  Vay.  +  vi  •  -1  -  Vai,  ■  “y,*i  !  J'i- 
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We  will  use  for  composition  of  equations  of  dynamics  of  an  air¬ 


craft  expressions  (3)  for  speeds  and  accelerations  of  elementary- 
points  of  the  aircraft  in  projections  on  connected  axes  of  coordinates. 
We  will  write  equations  of  forces  of  the  dynamics  of  the  system; 

n  it  /i  »  n  * 

(a).  *>,«,,  -V/>,  (6 ).  (a).  /m 

,_i  /-i  /-i  .  i-t  ' 


where 


are  projections  of  the  vector  of  ex¬ 
ternal  forces,  applied  to  the  system, 
n  is  the  quantity  of  elementary  particles 
of  the  aircraft  in  a  given  moment  of 
time . 


We  will  preliminarily  introduce  the  following  conceptions,  having 
a  real  physical  meaning,  which  we  will  use  in  the  future. 

1.  Speed  of  change  of  the  center  of  gravity  of  the  aircraft  at 
constant  mass. 

2.  Speed  of  change  of  the  center  of  gravity  of  the  aircraft  at 
variable  mass. 

3.  Speed  of  change  of  the  moment  of  Inertia  of  the  aircraft 
at  constant  mass. 

4.  Speed  of  change  of  the  moment  of  inertia  of  the  aircraft 
at  variable  mass. 


We  will  explain  this  examples; 

a)  let  us  have  on  the  aircraft  two  tanks  with  a  combustible 
liquid.  One  tank  -  in  the  nose,  and  the  other  -  in  the  tail,  and 
let  the  liquid  from  the  nose  tank  overflow  into  the  tail  tank.  The 
position  of  the  center  of  gravity  of  the  aircraft  (x  ,  y  )  changes 

eg  eg 

and  is  moved  to  the  tail.  The  Mass  of  the  aircraft  with  this  does 
not  change.  In  this  case  we  distinguish  the  speed  of  the  center  of 
gravity  of  the  aircraft  at  constant  mass  and  designate  correspondingly 


('«  —  const),  (in  —  const). 

The  moments  of  Inertia  of  the  aircraft  with  respect  to  axes  Xj,  ylt 
Zj  with  this  also  changes  and  the  speeds  of  their  change  will  be 
recorded: 

til,  v  ,  di. 

— -(in const),  —  ■  (///  const),  — -  (m  ^  const); 

ill  di  dt 

b)  let  the  combustible  liquid  continuously  overflow  from  the 
nose  tank  in  tall  tank  and  then  flow  from  the  aircraft  with  a 
determined  speed.  In  this  case  the  center  of  gravity  of  the  aircraft 
is  also  moved,  but  its  mass  with  this  is  variable.  In  this  case  we 
distinguish  the  speed  of  the  center  of  gravity  of  the  aircraft  at 
variable  mass  and  designate  as: 


tlx 

-2d  (//l 
dt 


var), 


(m  » 


var). 


The  moments  of  inertia  with  respect  to  axes  x,,  yi,  zlt  changing 
not  only  due  to  the  change  in  coordinates,  but  also  because  of  the 
change  of  the  mass  and  speed  if  their  change,  we  designate  corres¬ 
pondingly: 

,ii,  di,  di, 

- -tin  verb  -  '  (in  var),  - - (m  -=  var). 

,it  v  dt  dt 

In  the  future  these  physically  clear  conceptions  we  will  a1  so  use. 

For  brevity  of  spelling  of  derivatives  we  will  introduce  the 
following  conditional  designations: 


i 


!'  fill,  \ 

— -  (in  —  const)  —  (  — -  )  h  t.  a. 

dt  '  dt  •  mr 


Having  substituted  into  equation  (4a)  the  expression  for  w  , 

^  1 


we  obtains 


dV2 


dt 


j»i 

/-l 

l-J 

/i 

ft 

ft 

-  •><»,. 

» 

dy«  .  v\m  <***, 

u> 

Va:% 

Jtti,  -i-  u), 

Li 

*  dt  Li  dt* 

La  1 

<•1 

t- 1 

t  •  \ 

ft 

<• 

M 

II 

:1j 

«?.  VAyt  ltli  —  \ 

V 

»>,xu 

i  10  j.  V 

i-i 

i-1 

l—l 

/—l 

(5) 
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The  obtained  equation  (5)  can  be  simplified,  if  one  considers 
the  following  dependences: 

1.  The  sum  of  elementary  masses  of  the  aircraft  is  equal  to  the 
mass  of  the  entire  aircraft: 


V  /«,  /Vi 

l-l 


(  ba) 


2.  The  sum  of  the  static  moments  of  the  elementary  masses  of 
the  aircraft  is  equal  to  the  stati--  moment  of  its  whole  mass: 


V,fl.Z|/-=Mzll,.V,«(.y|,  Aly„„  £/«i-v„  *  Mx„ 

(_i  i-i  i~  i 


(6b) 


3.  The  sum  of  the  products  of  the  elementary  masses  of  the 
aircraft  on  corresponding  relative  velocities  is  equal  to  the  product 
of  the  entire  mass  on  the  speed  of  the  center  of  gravity  of  the 
aircraft  at  its  constancy: 


/c 


(6c) 


i-1 


/—I 


4.  The  sum  of  the  products  of  the  elementary  masses  on  the 
corresponding  relative  accelerations  is  equal  to  the  mass  of  the 
entire  aircraft  on  the  average  value  of  relative  acceleration: 

rf33'c~  ri  --  ( 6d) 

n,  — f-  =  m 
‘  ,ii' 
i-i 


=  M- X  =  M  -  Af  ^ . 

'  r/l>  dP  h  ‘  dP  dt> 


where  ZlfsL—  are  the  average  values  of  projections  of 

’  rfn  ’  dt- 

relative  acceleration  of  the  variable 
mass  of  the  aircraft  on  the  connected 
axes  of  coordinates. 

Let  us  note  that  at  a  constant  speed  of  supplying  of  liquid 
in  the  engine  the  average  values  of  projections  of  relative  accelera¬ 
tions  are  equal  to  zero.  Having  substituted  dependences  (6a),  (6b), 
(6c),  ( 6d)  into  equation  ( 5)  ,  we  obtain: 


/-I 


Analogous  to  (4a),  equations  (4b)  and  (4c)  take  the  following 

form: 


// 


1- 


M  ['->  +  •«.  •  - 2...  (%)  _  ,'jir ) 

l  ^  ^  ^  ^  'mf  K  lit  '  me 

+  rf/T-  +  +  “*'•  “'•  ’«  « ~  "i.-V,  T~  «x,  yn  t  -f  <»,, -t,  r  j  ■ 

-IN-  +  %  *  ,  -  ±£Xa  ,  +  2*»jr,  (<*■*)  _ 

‘  *tt  tit  \  rf<  /me 

A 

/rfjf,  B  T\  </;r 

~  2u’''  \~dr)„/'  'Ip  +  ,0t'  K*  ■*■  ",'  A'«  T  - ,o;,  *,  T  -  <■>.’■.  Vam,  - 

1  * 

—  WV.  u>y,^.ya  T  =£/a„ . 

J  /-! 


(8) 


We  obtained  equations  of  dynamic  equilibrium  of  the  variable 
mass  of  aircraft  in  projections  of  forces  on  axes  x1}  yi,  zj.  How¬ 
ever  these  equations  are  insufficient  for  description  of  the  conditions 
of  full  dynamic  equilibrium  of  the  system.  To  it  is  necessary  still 
to  add  equations  of  moments  around  axes  x1}  yi,  z,: 
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(9a) 

(9b) 

(9:) 


Having  substituted  into  equations  (9)  the  expressions  for  accel¬ 
erations  Wj,  w,,  W)  and  considering  with  this  that  the  mass  will  be 
displaced  symmetrically  relative  to  the  longitudinal  plane  of  the 
aircraft  and  that  the  speed  of  its  change  is  constant 
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we  obtain: 


{ttt  +  v4»  ~  }  +  ^  4« 1  + 

/-i 

/» 

2u,x'  “p  ~  (“x.  “»• - pp  /x^'.  +  wr.  w'.  (It,  ~  h)  — 

t-l 

v-l  /.I 

•  ("7T  +  Vm)  -  “A  ,  4-  ^  /,,  4-  + 

v-l 

rt 

+  2“.-,  pf* Zw  -  (l“v-  +  pr)  <»x.  (u  -it,)- 

/—l 

*  i  t  'J 


(  10a) 


(  10b) 


(~dT  ~  ^  _  t  +  —■!,,+  2<Uf,  yj"«i  ~  •*„+ 

K.  -  "»/.({, •.  -  /,,)  “(pjp*  +  “r.  -  “A  Va.*,  )  My„  T+ 

n  n 

(/r,v 


(  10c) 


Till  now  In  the  conclusion  of  the  equations  we  did  not  place 
limitations  on  the  selection  point  of  the  origin  of  the  coordinates 
of  the  connected  system.  In  the  selection  of  the  origin  of  coordi¬ 
nates  in  the  Immediate  center  of  gravity  of  the  aircraft  the  static 
moments  are  equal  to  zero  and,  consequently,  the  members  of  the  pre¬ 
ceding  equations  ( 7) ,  (8),  ( 10) ,  containing  static  moments  relative 
to  the  center  of  gravity,  vanish. 

In  the  preceding  equations  ( 10)  enter  also  members,  containing 


Vr,  ^'yu.  V««,  xu. 

'j  ■  Hi  *  1  ' ,  '  HI  “ 


These  sums  can  be  interpreted,  correspondingly,  as  derivatives  of  the 
product  of  inertia  according  to  time  at  a  separate  change  in  coordi¬ 
nates  x i ,  y i .  Having  differentiated  the  product  of  inertia  of  the 
entire  aircraft  according  to  time  at  constant  mass,  we  obtain: 

~t  It  1  at  ~L  ‘  <tt  iu  '  L  ‘  at  *u’  (  U' 

l~l  /-l  /-I 

We  introduce  the  following  designations: 

^  ( 12) 

2h  7r-v‘- =  7t  mc-  2/  -v»~  rfi  ",f- 

If  we  select  the  origin  of  the  connected  system  of  coordinates 
in  the  center  of  gravity,  considering  dependences  (  12)  and  (  13) 
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then  the  equations  of  motion  of  the  aircraft  ( 7)  ,  (8),  (  10)  take 
the  following  form: 
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In  the  right  sides  of  the  obtained  equations  are  contained  onl 
external  forces.  External  forces  and  moments,  acting  on  the  air¬ 
craft,  are: 

1)  aerodynamic  forces  and  moments  Xi,  Yi,  Zi,  M„  ,  ,  M  , 

x  1  / 1  z  i 

2)  gravity  G, 

3)  reactive  force  (  thrust)  $  and  moment  from  thrust  M  . 
Aerodynamic  moments  and  forces  at  linear  approximation  are 


expressed  by  the  following  dependences: 


MX,  “  ("C""'.  +  "CS-,  -■  s  -  /«*"«„  -f  mlfiSql, 
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These  coefficients  rn^,  m^x,  c  ,  c  ...  depend  on  M  numbers, 

y  i  x  i  y  x 

angles  of  attack,  and  also  on  the  position  of  the  center  of  gravity 
of  the  aircraft.  Since  the  equations  of  the  motion  of  the  aircraft 
are  composed  relative  to  the  system  of  coordinates,  passing  through 
the  center  of  gravity,  the  position  of  which  changes  in  time,  then 
it  is  natural  that  also  the  aerodynamic  coefficients  are  re-counted 
relative  to  the  current  position  of  the  center  of  gravity. 

Having  substituted  equality  (1?)  into  equations  ( 14)  ,  we  obtain 
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(16) 


The  obtained  equations  (6)  describe  the  motion  of  the  aircraft 

in  the  general  case  of  curvilinear  flight.  The  motion  of  the  aircraft 

is  interconnected  in  all  three  planes. 

During  longitudinal  motion  co  ,  cu  =  V  =  y  =  Zt  =  0  general 

y i  z  i 

equations  of  motion  (  16)  degenerate  into  the  equations  of  longitudinal 
motion  of  the  aircraft: 
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( 17) 


where 
gravity, 
and  the 


is  the  arm  of  reactive  force  relative  to  the  center  of 
;p  is  the  angle  between  the  direction  of  reactive  force 
axis  of  the  aircraft  (Fig.  4). 


Fig.  4.  Diagram  of  the  forces,  applied 
to  an  aircraft  in  the  longitudinal  plane. 
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If  one  takes  in  the  obtained  equations  ( l6) ,  ( 17)  the  relative 
velocity  of  fuel  mixture  of  the  aircraft  equal  to  zero,  then  they 
identically  coincide  with  the  equations  the  "solid  of"  aircraft. 
However  it  is  necessary  to  underline  that  although  the  obtained 
equations  remind  us  of  equations  of  motion  of  "solid"  aircraft, 
they  are  considerably  distinguished,  and  namely:  the  equations  of 
an  aircraft  with  a  liquid-fuel  rocket  engine  contain  mass  and  moments 
of  inertia  relative  to  the  connected  axes  -  variable  in  time,  the 
origin  of  coordinates  of  which  coincides  with  the  instantaneous 
center  of  gravity  of  the  aircraft  in  flight,  and  also  additional 
members,  considering  the  relative  motion  of  the  fuel  mixture  in  the 
aircraft. 

In  conclusion  we  will  spread  the  obtained  equations  of  dynamics 
of  an  aircraft  with  a  liquid-fuel  rocket  engine  in  the  case  of  a 
turbo-jet  engine. 

Analyzing  consecutively  our  computations,  it  is  easy  to  see  that 
we  nowhere  used  the  peculiarities  of  the  variable  mass,  concluded 
in  the  separation  of  combustible  gases  or  in  the  connection  of  air 
masses.  It  is  important  that  the  variable  mass  has  an  influence 
on  the  dynamics  of  the  aircraft  -  this  is  the  presence  of  relative 
velocity  of  particles  of  the  fuel  mixture  in  the  connected  system  of 
coordinates.  Consequently,  the  obtained  equations  are  completely 
spread  in  the  case  of  a  turbojet  engine  and  are  equations  of  an  air¬ 
craft  of  variable  mass  in  the  most  general  form.  However,  at  deter¬ 
mination  of  derivatives  of  the  moments  of  inertia,  the  coordinates  of 
the  center  of  gravity  of  the  aircraft,  entering  into  the  equations  of 
motion,  one  should  consider  the  peculiarity  of  that  or  any  other  jet 
engine,  i.e.  the  separated  and  joined  masses  <>f  the  aircraft.  Below 
are  listed  formulas  for  determination  of  the  <•  orrespondlng  derivative 

/S 


at  constant  mass  (  Para.  2)  . 

Analyzing  the  obtained  equations  (  16)  ,  ( 17) ,  it  is  possible 
to  come  to  following  qualitative  conclusions: 

a)  in  the  equations  of  forces  the  relative  velocity  of  the 
center  of  gravity  of  the  aircraft  functionally  have  influence  on 
the  force  of  inertia,  as  well  as  the  carrying  speed  of  the  center  of 
gravity, 

b)  in  the  equations  of  moments  of  speed  the  changes  of  the 
moments  of  inertia,  as,  for  example: 


have  an  influence  on  the  damping  qualities  of  the  aircraft.  The 

coefficients,  characterizing  their  own  damping  of  the  aircraft,  are 

m^*1,  m^1,  jn°x  1 ,  rnfz  1 ,  having  negative  numerical  values. 

*1  y i  y i  z i 

Consequently,  so  that  the  corresponding  speeds  of  change  of  moments 
of  inertia  they  improved  the  damping  qualities  of  the  aircraft,  it  is 
necessary  so  that  ,  f-^-\  would  be  positive  values, 

\  dt  /„V  V  >U  ),nc  \  <«  /  me 


and 


\  At  Jmc  \  At  J, 


negative. 


In  connection  with  these  conclusions  it  is  interesting  to  explain, 
under  what  conditions  the  corresponding  derivatives  of  moments  of 
inertia  have  positive  or  negative  values.  Before  being  occupied 
with  the  analysis  of  these  conditions,  we  will  introduce  formulas, 
determining  the  derivatives,  and  then  on  the  basis  of  the  formulas  we 
will  find  the  conditions  of  positiveness  and  negativity  of  these 
magnitudes. 

2.  Determination  of  speed:’.  oT  jhange  of  the  •enter_'jf  gravl  \:j 
and  moments  of  Inert. la  of  an  nlr'ran.  jit  ■•'umlaut  mass. 
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1.  Determination  of  speeds 
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me 

V  i it  ) 

\  <n  J 

In  the  determination  of  the  speed  of  the  center  of  gravity  of 
an  aircraft  (for  constant  mass)  we  will  originate  from  the  consider¬ 
ation  that  the  sum  of  static  moments  of  the  elementary  masses  of  an 
aircraft  is  equal  to  the  static  moment  of  the  entire  mass: 
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i-i 
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mtdy  Mdyn .. 
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The  left-hand  side  of  equation  ( l8a)  is  presented  in  the  following 

form: 


V  //!,(.*„  ,  t  x, 

Iwl 


n  n 

dx)  —  X  w,(.v,( ,  xx)  =  V  ,Uldx, 

t" I  /-I 


(19) 


where  xt  are  the  coordinates  of  masses  relative  to  the  system, 

passing  through  the  .'enter  of  gravity, 

X  is  the  coordinate  of  the  center  of  gravity  of  an  aircraft  rel- 
c6 

tive  the  original  system  of  coordinates. 

The  static  moment  of  an  aircraft  in  a  moment  of  time  t  +  dt 
relative  to  the  original  system  of  coordinates  is  equal  to 

m 

\m,{xa ,  +  •*!+■  dx).  w^ere  y„,  at  t  +  dt  are  equal  to  M2,  = 

(-i 

=  M2  +  dM.  In  order  to  obtain  the  static  moment  of  an  aircraft  in  a 
moment  of  time  t  +  dt  relative  tt  the  original  system  of  coordinates 
at  constant  mass,  i.e.  at  Mi,  it  is  necessary  to  add  to  the  last 
expression  the  static  moment  of  the  separated  mass  from  the  exhaust 
nozzle  of  the  engine  after  time  dt: 

V  m,  (xa  T  ~  .v.  i/x)  —  V  in,  (,v„  ,  x,  -i-  dx)  -f-  dMx  , 

“I'M 
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where  x„,  yu  are  the  coordinates  of  the  exhaust  nozzle  of  the  engine 
ri  n 

relative  to  the  original  axes  of  coordinates. 

Having  substituted  the  obtained  equality  in  ( 19)  ,  we  will  have: 

£  '«/(*«  T  4-  -V,  -f  tlx)  T  4-  x,)  -f  dMx%  =  Mdxa  , 
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or 
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Having  differentiated  in  the  preceding  equality  the  expression 
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,  we  obtain: 
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and  analogously 
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where  x0  =  x„  —  x  are  the  coordinates  of  the  exhaust  nozzle  of 

H  eg 

the  engine,  y0  =  yH  -  ycg  -  relative  to  the  systems  of  coordinates, 
passing  through  the  center  of  gravity. 

Spreading  the  obtained  equality  on  the  aircraft  with  a  turbo¬ 
jet  engine  (TRD) ,  when  besides  the  division  ;f  masses  takes  place 
also  their  connection,  we  obtain: 
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dt  v  d: 
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where  __  is  the  second  expenditure  of  separated  masses, 

dt 

dm 

•it  is  the  second  expenditure  of  joined  masses. 

If  one  considers  that  for  an  aircraft  with  a  turbojet  engine  with  an 
accuracy  sufficient  for  practice  it  is  possible  to  take 


T 

dt 


( m  =  var) 


( m  -  var)  ~  0, 


dm 

_  rf.U 

dt 

' dt' 

JL/ 


then  we  find 
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We  obtained  simple  formulas  for  determination  of  the  speed  of  the 
center  of  gravity  of  an  aircraft  at  m  =  const. 

2.  Determination  of  the  speed  of  change  of  moments  of  inertia 
of  an  aircraft  relative  to  the  center  of  gravity. 

Let  during  the  time  dt  the  coordinates  of  masses  mi  relative 
to  the  system  of  coordinates,  passing  through  the  center  of  gravity, 
be  changed  to  dx,  dy.  We  find  expressions  of  moments  of  inertia  of 
the  aircraft  with  respect  to  the  axes,  passing  through  the  center 
of  gravity  at  constant  mass.  The  moment  of  inertia  of  the  aircraft 
I  in  a  moment  of  time  t  +  dt  (with  a  decreased  mass  on  dM)  relative 

x  1 

to  the  former  position  of  the  center  of  gravity  is  determined  by 
the  formula: 


.\Uy\ 


where  mass  M2  corresponds  to  time  t  +  dt.  So  that  we  obtain  a  moment 

of  inertia  of  the  aircraft  in  moment  t  +  dt  relative  to  the  former 

center  of  gravity  with  mass  Mi,  it  is  necessary  to  value  I  ( M2)  to 

•X  1 

add  its  own  moment  of  inertia  of  the  discarded  mass  of  fuel  and 

magnitude,  equal  to  the  product  of  the  second  expenditure  of  mass  of 

fuel  on  the  square  of  coordinates  of  the  output  section  of  the  engine 

relative  to  the  center  of  gravity  and  on  dt. 

Thus,  the  moment  of  inertia  I  of  mass  Mi  at  t  +  dt  is  determined 

x  i 

by  the  following  formula: 
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where  r^  is  the  radius  of  gyration  of  a  section  of  the  discarded 
part  of  *-he  fuel  mixture. 

To  find  Increase  of  the  moment  of  inertia  I  of  mass  Mi: 

x  i 

4.  dt)  -  / .  (/)  =  /liU,  ,,rj  -  M,dy\ ,  -  tIMyl  +  </AIr;  -  /r< „  (AI,). 

Having  divided  the  preceding  equality  by  dt,  we  obtain: 
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As  can  be  seen  from  the  formula,  in  order  to  determine  the 
derivative  of  the  moment  of  inertia  in  time,  it  is  necessary  from  a 

graph  I  =.  f  ( t)  to  determine  the  derivative  and  to  add  to  it  the 
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where  x„,  y„  are  the  coordinates  of  the  exhaust  nozzle  of  the  engine 

n  H 

relative  to  the  original  axes  of  coordinates. 

For  aircraft,  supplied  with  turbo-jet  engines  ( TP.D)  ,  in  a  case, 
when  takes  place  division  of  combustible  gases  and  connection  of 
air  masses,  the  preceding  equalities  take  the  form: 
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(m  =  const)  --  j  "L \(V!  x  r*  i  _ \i!L  I  {f-  x  f3 ),  4  . 

•It  V  ,11  ) nV  !  dt  !-"  r  |  dt  I  ^  '  *  *m* 

Y  ’  <-  -  -  Cl’  !).,+!  «  'J  >« -Ifi «  *  f  W 

■(-r-v.Y'fY  w%,~ 

—  ■  j.  r1\ 


£3 


where  _  is  the  second  expenditure  of  separated  mass, 

it 

dm 

is  the  second  expenditure  of  joined  mass. 


If  one  considers  that  for  an  aircraft  with  a  turbojet  engine 
exist  the  dependences 
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then  the  preceding  equalities  can  be  presented  as: 

(m  ~  const)  -  |  1  |(jr2  +  r;  ),*,--  ( y\  +  r\  )(ln)). 

(„  =  const)  -1^1 1(4  ~  r%,-  (*l  +  r*  )(m,|. 

— (/»  —  const)  --  j  [{.<1  -  >’5  -f  rpIM)  -  (*ij  +  vj  +  r;)(„)|. 
it  ,  it  \ 


The  obtained  formulas  are  very  convenient  to  use  for  determination 
of  the  speeds  of  moments  of  inertia  of  an  aircraft  at  constant  mass. 


necessary  for  calculation  of  the  dynamics  of  the  aircraft. 


Determination  of  (d,*.y 


/it x,  v,  in  \ 

V  it  J, 


f  (y)\ 

me  ^  * 


For  determination  of  derivatives  of  the  product  of  inertia  in 
time  at  a  separate  change  of  coordinates  it  is  necessary  to  directly 


calculate  in  the  function  of  time  any  of  the  two  partial  derivatives. 

.  itx  v,  (•«)\  i'  ,  l ix 

Determining,  for  Instance  ■ — — j  •=  J  tlnl~t  J/>  ,  the  second 

mt  (m) 


partial  derivative,  we  obtain  by  means  of  subtraction  from  the  total 


derivative 


t'Hx,  vMh 
\  it  J  me 


We  will  show  a  method  of  determination  of  the  total  derivative 


of  the  product  of  inertia  at  constant  mass.  In  the  moment  of  time 

t  +  dt  the  product  of  inertia  is  equal  to  I  ( M2)  ,  and  relative  to 

x  iy  i 

the  system  of  coordinates,  passing  through  the  center  of  gravity  in 


the  moment  of  time  t,  it  equals  I  ( M2)  +  M2dx,  •  dy 

iy  i  eg 

The  product  of  inertia  in  moment  of  time  t  +  dt  from  mass  Mt 
equals  the  sum  of  the  product  of  inertia  from  mass  M2  and  the  product 
of  the  separated  mass  after  time  dt  on  coordinates  x0,  yo  of  the 
exhaust  nozzle 

/>uVl(AU  =  /„,(•«,)  +  ,  +  dMx,y0. 

We  will  record  the  Increase  of  the  moment  of  inertia  after  time  dt: 
dI*sSm  ”  const)  "  Mt<ixa,dyu  ,  -  tIMxay0. 


Dividing  the  increase  by  dt,  we  have: 


rf/ 


-(>«  =  const) 


u  t  \ 

rii  A 


J.U 

dr 


xoy<>- 


(26) 


For  an  aircraft  with  a  turbojet  engine  the  preceding  equality 
takes  the  form: 


(//i 


const)  “  ^ 


If  one  considers  that 


,ur 


dt 


!  <A,v,  u  t 

1 M  i  v 

\  dt 

/'in  r 

!  dt  b 

/A, 

3  dt" 

0  - 
J  mV 

0,  \^  \ 
1  dt  | 

const)  =*  j 

(Ml1, 

Til' 

•v')j'o(‘'i)  —  j -/,n 

I  (/<  ; 


(Ml 


Thus,  in  order  to  determine  the  total  derivative  from  the 

product  of  inertia  in  time  at  m  =  const,  it  is  necessary  to  determine 

from  a  graph  of  the  product  of  inertia  relative  to  the  system  of 

coordinates,  passing  through  the  center  of  gravity,  the  time  derivatlv 

and  to  add  the  product  \~- 

I  dt 

3 •  Evaluation  of  the  Influence  of  variability  of  mass  In 
the  motion  equations  of  an  aircraft. 


Equations  of  mo!  Ion  of  an  aircraft  of  variable  mam  are 


distinguished  from  analogous  equations  of  motion  of  a  solid  aircraft 
into  every  moment  of  time  by  additional  members,  considering  the 
speed  of  change  of  the  center  of  gravity  and  the  speed  of  change 
of  the  moments  of  inertia  of  the  aircraft  at  constant  mass.  Above, 
on  the  basis  of  analysis  of  equations  of  motion  of  an  aircraft  (  l6) , 

we  arrived  to  the  conclusion  that  for  improvement  of  damping  it  is 

\  fd,y>  \  would  be  positive 

necessary,  in  order  to  {  ,u  )J  \  At  jmc 

values. 

Ho»  we  will  define  the  qualitative  conditions,  with  which  these 
requirements  are  executed.  For  this  purpose  we  will  turn  to  the 
formulas,  determining  the  speed  of  change  of  the  moments  of  Inertia 

at  constant  mass. 

Let  us  consider,  for  example,  the  speed  of  change  of  the  moment 
of  inertia  around  axis  yx,  determined  by  the  formula: 


!dl±_\  _  (Al±\  r\). 

\  dt  )„r  \  dt 


As  can  be  seen  from  the  last  formula  so  that  )*').,  iS  SrL'ator  than 
zero,  fulfillment  of  the  following  condition  is  necessary: 


!£!  I  ( jl  r? )  >  — 

it  |l  0  •  dt 


dt  \mV 


Fulfillment  of  this 


condition  is  possible  in  two  versions. 


w 


Fig.  5.  Diagram  of  the  location  of  the 
engine  in  front  of  the  center  of  gravity 
of  an  aircraft.  l)  Fuel  mixture;  2)  center 
of  gravity. 


First  version:  the  exhaust  nozzle  of  the  engine  is  located 
In  front  of  the  center  of  gravity  of  the  aircraft,  and  the  center 
of  gravity  of  the  fuel  mixture  is  located  relative  to  the  center  of 
gravity  of  the  aircraft  at  significantly  smaller  distance  (Fig.  5)  • 


Fig.  6.  Diagram  of  the  location  of  the 
engine  behind  the  center  of  gravity  of  an 
aircraft.  l)  Fuel  mixture;  2)  center  of 
gravity. 


Second  verson;  the  exhaust  nozzle  of  the  engine  is  located 
behind  the  center  of  gravity  of  the  aircraft  at  a  large  distance, 
and  the  center  of  gravity  of  the  fuel  mixture  relative  to  the  center 
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of  gravity  of  the  aircraft  at  a  large  distance,  and  the  center  of 
gravity  of  the  fuel  mixture  relative  to  the  center  of  gravity  of 
the  aircraft  is  located  at  a  significantly  smaller  distance  (Pig.  6). 

-0.0} 

use 

cot 

••  «  4  4  -V  '•> 


Fig.  7. 

In  the  arrangement  of  an  aircraft  one  whould  consider  these 

/d /.,  \ 

considerations.  We  will  note  however,  so  that  >0. 

is  actually  a  realizable  arrangement  according  to  the  second  verson. 
Analogous  reasonings  are  spread  also  on  the  speed  of  change  of  the 
moments  of  inertia  around  axes  Xj,  zx.  The  above-introduced  analysis 
has  a  qualitative  character.  However  in  practice  it  can  appear  that 
the  calculation  of  additional  members  numerically  have  little  Influence 
on  equations  of  motion  of  aircraft.  Therefore  it  is  expedient  to 
produce  a  numerical  evaluation  of  additional  members  for  a  judgement 
about  the  boundaries  of  application  of  these  equations.  Let  us 
consider  an  aircraft,  accomplishing  takeoff  with  the  help  of  a 
liquid-fuel  rocket  engine  from  the  method  of  horizontal  flight. 
Aerodynamic  and  flight  characteristics  of  an  aircraft  during  take 
off,  as,  for  example,  "C'\  •  "h,!  •$<//,  Sqba,  arc  shown  in 

graphs  ( Fig.  7)  • 

The  inertial  weight  characteristics  of  an  aircraft  during  take 


;?/ 


off  In  a  function 


of  time  are  given  on  the  gt’ 


*aphs  ( Fig-  8)  . 


Differentiating  In  time  the  curves  of  change  0,  V  H,’  X=s’ 

v  „e  wlll  find  the  speed  of  change  of  these  magnitudes  at 
yc%’ 

variable  mass,  and  namely: 

,i  /til,  \  ilM  S'11'.',' 

&V  &)-•  * ' ( '** 

Using  then  formulas  (25).  (25).  <  *>  •  «detemlno  the  speeds  of 
Lge  or  the  center  of  gravity  and  moments  of  Inertia  of  an  aircraft 

at  constant  mass: 


\  At  9K 


•  *j\ 

1  It 


\  tit  /v 


w.,L'  v « J 


On  the  basis  of  data  of  the  numerical  calculation  of  the  above- 


considered  aircraft  we  will  produce  calculation  of  coefficients  of 
additional  damping  from  the  variability  of  mass  of  the  aircraft. 

Prom  equations  of  motion  of  an  aircraft  ( l6)  it  follows  that 
coefficients  of  additional  damping  of  an  aircraft  from  variability 
of  mass  are  determined  by  the  following  formulas: 


-  f  —  - 
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A/// 
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Results  of  the  calculations  of  these  coefficients  are  given 
in  the  table. 
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In  comparing  the  numerical  values  of  additional  coefficients 
of  damping  with  its  own  coefficients  of  damping,  we  notice  that  in 
the  equations  of  moments  ( l5)  the  additional  members 


can  considerably  have  an  Influence 


(d,*>  \  \ 

V  <tt  )mt’  \  it  Jm.’  \  ,it  Jmc  '  V  it'  ), 

on  the  damping  qualities  of  the  aircraft  and  it  is  impossible  to 
disregard  them.  In  projection  of  an  aircraft  one  should  consider 
this  circumstance.  We  will  note  however,  that  in  the  equations  of 
motion  of  an  aircraft  one  whould  consider  the  coefficients  of  addi¬ 
tional  damping,  when  they  exceed  10$  in  reference  to  its  own  damping 
of  the  aircraft,  since  the  aerodynamic  coefficients  are  usually 
given  with  an  accuracy  to  10$. 

In  the  equations  of  forces  ( l6)  it  is  possible  to  disregard  the 
members,  considering  the  speed  of  change  of  the  center  of  gravity  at 
constant  mass  ( ^JLl\  ,  (— - \  ,  in  view  of  their  insignificant 

\  J  me  \  <it  Jmc 

influence. 

In  conclusion  we  will  consider  the  case  of  flight  of  an  aircraft, 
when  its  mass  changes  abruptly,  for  example  due  to  booster 
separation.  Do  there  appear  in  this  case  additional  damping  coef¬ 
ficients  in  the  equations  of  motion?  We  will  show  that  the  speeds 
of  change  of  moments  of  inertia  at  constant  mass  in  this  case  equals 
zero.  Let  us  consider,  for  example,  the  speed  of  change  of  the 
moment  of  inertia  of  an  aircraft  with  respect  to  axis  y i  —  \~%~j 
at  booster  separation. 

Let  its  own  moment  of  inertia  of  the  booster  equal  m^r®, 
where  myc  is  the  mass  of  the  booster, 
ry  is  the  radius  of  gyration. 

Then  the  speed  of  change  of  the  moment  of  inertia  of  the  aircraft 
relative  to  the  center  of  gravity  of  the  aricraft  at  variable  ma33 
will  be: 


\  it  jr.V 


”yc  Uj  T  ry  ) 

it 
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The  speed  of  change  of  the  moment  of  inertia  of  the  aircraft  at 
constant  mass  is  determined,  as  it  is  known,  by  the  following  depen¬ 
dence  t 


(dM 

V  dl  V  dt  JrnV  « 


(A/ 


0). 


Substituting  the  preceding  in  the  obtained  equality,  we  will  have: 

(d,y . '  ^  ~ '"■>< ^  -  '•y)  4-  %(*'!-'  4 >  _  o)  =  0 

V  dt  )mc 

Consequently,  we  showed  that  the  speed  of  change  of  the  moments 
of  inertia  in  this  case  is  equal  to  zero  and  therefore  in  the  equa¬ 
tions  of  motion  members  of  additional  damping  are  absent.  However, 
the  value  of  mass  should  be  new  and  the  aerodynamic  coefficients 
calculated  correspondingly  to  the  new  position  of  the  center  of  grav¬ 
ity.  Equations  of  motion  of  aircraft  are  solved  under  the  original 
conditions,  corresponding  to  the  moment  of  booster  separation. 

-  '  c„  .  <  ,  r  / 

Conclusions 

On  the  basis  of  what  has  been  stated  it  follows  that  the  motion 
of  an  aircraft  with  a  jet  engine,  thd  mass  and  moments  of  Inertia. of 
which  greatly  vary,  should  be  depicted  with  equations  ( 16)  ,  considering 
the  variability  of  mass.  If  the  mass  is  variable  and  with  this  takes 
place  the  relative  motion  of  the  fuel  mixture  in  a  connected  system 
of  coordinates  of  the  aircraft,  then  in  the  equations  of  motion  of 
the  aircraft  is  considered  the  variability  of  the  mass  and  the  addi¬ 
tional  members,  considering  the  damping  influence  on  aircraft  of 
relative  motion  of  fuel  mixture. 

It  is  very  important  in  the  projection  of  an  aircraft  to 

C 

correctly  arrange  the  position  of  the  engine  (booster)  and  tanks  of 
fuel  mixture  .relative  to  the  center  of  gravity  of  the  aircraft.  For 

?73-yT-6i-9L  /U-2  ££ 


yf*L  *  c.o- 

an  increase  of  damping  of  an  aircraft  it  is  expedient  to  -dispose  the 

Jet  engine  behind  the  center  of  gravity  at  a  significant  distance. 

Additional  damping  from  variability  of  mass  one  should  consider  in  'a 

case,  when  the  coefficients  (—)  — — 155_  etc. 

V  dt  Jmc  Sql  m».r,  \  dt  m-r, 

are  greater  than  10^.  *' 

In  the  equations  of  forces  ( l6)  it  is  possible  to  disregard  the 
members,  considering  the  relative  motion  of  the  center  of  gravity  of 
an  aircrafts 


(dx"\ 

(d*"\ 

V  dt  )me' 

\  dt  JJ 

If,  however,  the  mass  of  an  aircraft  is  variable,  but  with  this 
there  is  no  relative  motion  of  separated  masses  in  the  connected 
system  of  coordinates  (booster  separation) ,  then  the  members  of 
equations  (  l6)  ,  considering  the  additional  damping,  equal  zero.  In 
this  case  equations  ( l6)  consider  only  the  absolute  change  of  mass, 
the  moments  of  inertia  and  the  change  of  the  center  of  gravity  of  the 
aircraft. 
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